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We propose a new relation between quark mixing matrix and lepton mixing matrix. 
Since the parameters in the quark sector are well determined, we employ them to de- 
scribe the mixing of leptons. Phenomenologically, we study the neutrino oscillation 
probabilities for different channels, which can be measured precisely in forthcom- 
ing reactor and accelerator experiments. As an example of the applicability of our 
assumption, CP violation in the lepton sector is also discussed. In the latest T2K 
experiment, the range of the mixing angle ^13 is measured, and our prediction of 
^13 is compatible with their result. 
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The results of many neutrino oscillation experiments in the last decade have 
demonstrated that there exist physics beyond standard model in neutrino 
sector. It is commonly accepted that neutrinos are massive and mixing [1]. 
Neutrino oscillation is governed by two mass square differences Am2i, Am|]^ 
and the lepton mixing matrix proposed by Pontecorvo, Maki, Nakawaga and 
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If neutrinos are of the Majorana type, there should be an additional diago- 
nal matrix with two Majorana phases P = diag(e*"^/^, e*"^/^, 1) multiplied to 
Eq. (II]) from the right. But the two Majorana phases do not affect neutrino 
oscillations, thus we do not include them in our calculations. 

Before more underlying theory of the origin of the mixing is found, parametriz- 
ing the mixing matrix properly is helpful in understanding the mixing pattern 
and analyzing experimental results. A commonly used form is the standard 
parametrization proposed by Chau and Keung (CK) [3] 
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where s 



1-3 



sin^jj and 



cos 6 ij = 1,2,3) are the mixing angles, and 
S is the CP-violating phase. The current global fits of the neutrino mixing 
parameters are given at the l(3)cr level by [3] 
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These results for angles are compatible with the tri-bimaximal (TB) matrix [S] 
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Therefore, it is widely accepted that U^b is a good approximation to reality [6]. 

In contrast to the large mixing in lepton sector, the observed Cabibbo-Kobayashi- 
Maskawa [7](CKM) quark mixing matrix Vckm is close to the unit matrix. 
Although it seems that the mixing of quarks and leptons are unrelated with 
each other, there indeed exist phenomenological relations between mixing an- 
gles called quark-lepton complementarity (QLC) [Sj, given by 



e?2 + O12 = 7r/4, eg + 023 = vr/4, ~ 6 
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We pointed out in previous works ^ that the unit matrix pattern for quark 
mixing is connected with the bimaximal matrix pattern for lepton mixing 
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through QLC relations. But present data imply that the tri-bimaximal matrix 
f/xB is closer to reality than bimaximal matrix. Therefore, a natural idea is 
to connect the unit matrix in quark sector with the tri-bimaximal matrix in 
lepton sector. Based on this consideration, we propose here a simple relation 
between the lepton and quark mixing matrices 
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under an Ansatz that t/pMNS becomes an exact tri-bimaximal mixing U^-q 
in a limit Vckm = 1 [H]- Such a relation maybe comes from certain flavor 
symmetries, and the corresponding symmetry breaking effects may induce 
the deviations of the observed Vckm and f/pMNS from the exact unit matrix 
and tri-bimaximal matrix. By assuming Eq. ([7]), we can employ one set of 
parameters to describe both quark and lepton mixing matrices, thus Eq. ([71) 
can be regarded as the quark-lepton complementarity in matrix form. 



Currently, the quark mixing matrix is well determined. For the Wolfenstein 
parametrization 



CKM 



X 



-X 1 



AX^{p — if]) 
AX^ 
1 



0{X' 



(8) 



the up-to-date fit gives [13 



A = 0.2253 ±0.0007 , 
p(1-AV2 + ...) = 0.132; 
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By inserting Eq. ([8]) into Eq. one can easily get the PMNS matrix in terms 
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of A, A, p and 77, given bjfll 
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Clearly, we describe the deviations of Vckm from the unit matrix and ?7pmns 
from the tri-bimaximal matrix with the same set of parameters. In other words, 
relation Eq. ([7]) provides a unified way of parametrizing both quark and lepton 
mixing matrices. 

The lepton mixing angles are given by 
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With the central values in Eq. (jH]), we get numerically 



012 = 31.80°, = 44.66°, 01.3 = 7.67 



^23 
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which are compatible with the fit results in Eq. ([3]). Comparing with the exact 
tri-bimaximal mixing angles 



012 = arcsin 023 = vr/4, 0i3 = 0, 



(13) 



our result for 0i2 is more closer to the prediction of QLC since the la range 
for the corresponding quark mixing ang le reads [14J 0i^2 = 13.03(±0.05)°. 



For simplicity, we present expressions to the second order of A here, but all the 
results below result from C/pMNS to C'(A^). 
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Denoting e^j the deviations from the exact tri-bimaximal angles, we have 

ei2 = 012 - arcsin = -3.46°, £23 = 623-^ = -0.34°, €13 = ^13 = 7.67°. 

(14) 

As it shows, a relatively large 6*13 is predicted from our assumption, and such 
a prediction can be tested precisely in future reactor and accelerator experi- 
ments. 

An important property of the relation in Eq. ([7]) is that different phase con- 
ventions of VcKM give different predictions on lepton mixing angles. We give 
a brief argument here. If Vckm is rephased by taking 

^CKM = ^l^CKM^L (15) 

in which = diag(e*", e"^, e**) and "^2 = diag(e*'^, e**, e**) consist of arbitrary 
phases and can be absorbed into the redefinition of quark phases, Eq. ([7j) turns 
into 

Because Vckm (^ckm) ^'^^^ commute with ^2 (^2), it is generally impos- 
sible to absorb the phases in \l/2 (^2) i^^to the redefinition of lepton phases. 
As a result, the magnitudes of elements in ?7pmns and consequently lepton 
mixing angles depend on the rephasing matrix \l/2 (^1^2)) which does not bring 
any difference in Vckm- Therefore, predictions on lepton mixing resulting from 
Eq. ([7]) will be changed if we change the parametrization of quark mixing ma- 
trix. Similar arguments can be applied for assumptions discussed in Ref. [TT] . 
For generality, detailed discussion concerning the behavior of mixing angles 
on phases in \l/2 (^2) is given in the Appendix of this letter. As an instance, 
another parametrization [15j for Vckm is also employed to determine leptonic 
mixing angles in the Appendix. 

We now turn to the application of the assumption Eq. ((Tj) to neutrino oscil- 
lations. Let us denote Pq,/3 = Pii^a — ^ '^/s) the probability of transition from a 
neutrino flavor a to a neutrino flavor /3. Similar to the discussion in Ref. [T6] . 
the probability can be found as Pap = |5'/3a(t, to)p5 in which S{t,to) is the 
evolution matrix such that 

\u{t)) = S{t,toMto)), S{to,to) = l. (17) 

For simplicity, we neglect the effects due to interactions between neutrinos and 
matters in which the neutrino beam propagates and only deal with oscillation 



5 



probabilities in vacuum, the evolution matrix can be given by 



3 

S^ait, to) = T.(Ua^yUf,^e-''''\ a, (3 = 6, /i, T, (18) 

1=1 

where L = t — to is the length of the baseline in neutrino experiment and Ei 
are the eigenvalue of the effective hamiltonian 

H ~ ^f/diag(0, Ami,, ^^li)U\ (19) 
in which E is the average energy of the neutrino beam. 

Inspecting the values of the mass square difference in Eq. (jlj), one can find 
Am^i is much less than Am'^i, thus can be neglected to a good precision. The 
calculation of Pa/3 is now straightforward, such that by combining Eq. ([7]), 
Eq. dHD and Eq. (ITS]) , we get the evolution matrix S'^q, and consequently the 
oscillation probability Pap- Expanded to A^, oscillation probabilities can be 
expressed in a matrix form as 



P 



1 - 2A2A + 4yf AA^A X'^A-2^AX^A A^A - 2yf AA^A ^ 
A^A - 2 AA=^A 1-A A - A^A + 2 AA^A 



A^A - 2 JfAA^A A - A^A + 2JlA\^A 1 - A 



(20) 



where we have defined A = sin^ (-^^^ 



4_E 



Let us have a first look at the oscillation probability matrix in Eq. fl20|) . Appar- 
ently, it exhibits a hierarchical structure among different channels of neutrino 
oscillation. The diagonal elements Pee, P^i^i and P^r which measure the dis- 
appearance probabilities for z/g, and z/^ beams are of 0{1). It is also not 
difficult to understand that the <H- Vj. probabilities P^^ and Pr^ are of 0{1), 
since and are maximally mixing implied by data. Other terms measuring 
probabilities of Ue ^ z/^ and Ve ^ are of C(A^). Interestingly, there are no 
terms proportional to A, that is because such terms are suppressed by Am\,, 
which we neglect here. 

To get the anti-neutrino oscillation probabihties P^^, one needs to take the 
replacement t/pMNS ~^ ^pmns; which, in our case, is equivalent to reverse the 
sign of r], i.e. 

P-a-p = Po.fiij]^-r\). (21) 
However, the CP- violating parameter i] is missing in Eq. (!20ll . meaning that 
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to this approximation, we have 



al3 



(22) 



Therefore, CP symmetry is preserved, and resulting from the CPT theo- 
rem, time reversal symmetry is also preserved, i.e., Pap = Ppa- Consistently, 
Eq. f l2Up is a symmetric matrix. 



Defining the asymmetries in neutrino oscillations as 
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CP 



PnR — P; 
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(23) 



the probabilities in Eq. fl20p imply = 0. The reason for the vanishing of 
A^J is that we neglect the contribution of the smaller mass square difference 
Am^i. Now taking this into account, we get nonzero values for CP asymmetries 

as 
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in which 

a = -(2v^-3)Ar/AMsin(^^ 



sm I 



,L{Am 



Ami,] 



2E 



sm 



2E 



(25) 



The structure of asymmetries in Eq. comes from the unitarity of PMNS 
matrix and the conservation of probability. One also has [1] 
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in which 



Jcp = im{u^,Ue2u:,u; 



fi2) 



(27) 



is the rephasing invariant [18]. Combining Eq. f l2^ . Eq. (|25|) and Eq. (I26l) . we 

get 



■^CP = -i(2v^ - 3)Ar/A3 + O(A^). 

io 



(2J 



Compared with the results in Ref. |TT], where they arrive at Jcp ~ by 
assuming VbKMt^PMNS = Utb or V^^y^Upum = Utb, prediction of Jcp here 
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is quite smaller. This is because that the exact Utb implies Jcp = 0, thus 
Jcp only depends on the deviation of ?7pmns from [/tb (which we denote by 
D below). Since Vckm is close to the unit matrix, it can be regarded as the 
measurement of D. Then it is not difficult to understand the smallness of Jcp 
in Eq. as one has 

Dr^OiVcKu) and Jcp ~ C(A) (29) 

in Ref. [II], while 

^~(^(^ckm) and Jcp~0(A=^) (30) 

in this letter. 

We emphasize that in the very recent T2K experiment ^19j, observations of 
i^fj. oscillation indicate that at 90% C.L., the data are consistent with 

0.03(0.04) < sin^ 2^13 < 0.28(0.34) (31) 

for 6 = and normal (inverted) hierarchy. Such a result implies an apparent 
deviation from Utb and is important in the establishment of lepton mixing 
pattern, in which new symmetries among leptons may hide. If ^^13 is really 
large, the test of CP violation in the lepton sector is possible for future neutrino 
experiments since 6 is always multiplied by 6^13 in t/pMNS- With Eq. (!3T]) . 
straightforward calculations give 

4.99°(5.77°) < 013 < 15.97°(17.83°) (32) 

at 90% C.L., for (5 = and normal (inverted) hierarchy. We find that our 
prediction, i.e., ^13 = 7.67 is compatible with the T2K result. Thus the relation 
Eq. ([7]) may serve as a good description of the deviation of t/pMNS from Utb- 

In summary, we propose a new relation between quark mixing matrix Vckm 
and lepton mixing matrix f/pMNS given by Eq. ([7]) as the quark-lepton com- 
plementarity in matrix form. Based on this relation, we parametrize ?7pmns 
with the quark mixing parameters in Eq. f fTOl) . and determine the deviations 
of mixing angles from the exact tri-bimaximal angles. Especially, our predic- 
tion of the mixing angle ^13 agrees with the latest T2K result. For neutrino 
oscillations, we derive oscillation probabilities for different channels given by 
Eq. f l20|) . As we can see, to a good precision, the expressions for each proba- 
bility are quite simple. Furthermore, the CP violation in neutrino flavor tran- 
sitions are discussed. 

This work is partially supported by National Natural Science Foundation of 
China (Grants No. 11021092, No. 10975003, No. 11035003) and by the Key 
Grant Project of Chinese Ministry of Education (Grant No. 305001). 
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Appendix : Dependence of lepton mixing angles on the phase con- 
vention of CKM matrix 



With discussions concerning Eq. f|T6|) . we liave pointed out that ttie predicted 
lepton mixing angles are dependent on the phase convention of Vckm, i-e., 
\l/2 (^1/2) in our notation. By substituting the Wolfenstein matrix Eq. ([8]) into 
Eq. f|T6|) . we arrive at the general expression of t/pMNS, which includes explicitly 
the phases in and \l/2- Since phases in can be absorbed into the 

redefinition of lepton fields, the physical mixing angles are thus dependent 
only on phase parameters d, s and b. To a good accuracy, mixing angles are 
determined by 



sin^^ 
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sm U23 



12 1 

- - -A fv^cosa -cosa) + -A^ (2^2 cos 2a - 4^2 + s) + O(A^), 



sin^^ 



13 
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^ + A2 ^-^Acos(3 + Acos(3-^ \ +0(A3^ 
y-yf^A=^cos7 + 0(A^) 



where we define a = d~s, (3 = b — s and 7 = b + d — 2s, and the last 
equation is to 0{\^) because of the smallness of 613. It is easy to see that 
the sensitivities to phases of sin^6'i2, sin^ 6*23 and sin^ 6*13 are of C(A), C(A^) 
and C(A^) respectively. With the best fit values in Eq. ([9]), the dependence 
of each mixing angle is illustrated in Fig. [Hand Fig. [21 which show that, the 
prediction of 6*12 is strongly dependent on the phase, such that some areas of 
a are excluded by data and, the result in Eq. ( [T2l) can be improved by choose 
a particular value for a. However, the dependence of the predicted ^23 and ^13 
on phases is negligibly small compared with current data. 

Sin^ 623 




0.60 

0.55 



0.45 

0.40 
0.35 



Fig. 1. Behavior of sin^ 612 versus a = d—s (left) and sin^ ^23 versus (3 = b—s (right), 
both from — tt to tt. The soUd horizontal lines denote the So" ranges calculated from 
Eq. ([3|), while the dashed lines denote the best fit values. 



In order to demonstrate that our predictions of f/pMNS and mixing angles 
depend on the parametrization of Vckm, we here employ another Wolfenstein- 
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Sin2ei3 



0.06 



-3 -2 -1 

Fig. 2. Behavior of sin^ ^13 versus 7 = 6 + d — 2s from — vr to vr. The sohd (dashed) 
horizontal Hnes denote 90% C.L. ranges of the newest T2K result with normal 
(invertd) hierarchy. 

like parametrization [13], given by 
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to deduce the form of f/pMNS- The ranges for the parameters are 



A = 0.2253 ± 0.0007, h = 0.303t[!:|Ji^, 



/ = 0.754^°:°?2^ 5« = 90.97°+'-''' 



4.440 . 



Further detailed analysis and discussions concerning the relationship of this 
form of parametrization with others are given in Ref. |20]. By substituting 
this new Wolfenstein-like parametrization into Eq. (^^, we get the expression 
for [/pMNS to C(A^) as 



f^PMNS = yCKuUTBycKM = UtB + ^ 
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Comparing with Eq. ( fTOl) . one can easily find that differences begin to appear 
in terms of (9(A^). To 0{X^), mixing angles are given by 
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e e'" ; 

^/3 / 



sin^ 6*13 = ^ - \/ o (/ + COS 5) A' 



sm 



^23 = ^ - 3^ (4(^6 - 3)/ + 4(V6 - 3h) cos <5 + 3) A^ + ^ ((^2 - 2)/ + ^2 cos A^, 



sin2 6' 



12 



2\/3 

/2 - 1)A + ;^(3 - 2V2)A2 + i (2^/6/ + (3(\/6 - 2)/i - VQ) cos 5 + 9\/2 - I2) A^ 
3 3 D y V / 
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in which sin^ ^13, sin^ 9i2 differ from Eq. (fTTj) only in terms proportional to 
and the difference of sin^ ^23 is of 0{X^). Therefore, the numerical results for 
mixing angles are very close to Eq. f|T2|) . as straightforward calculation gives 

9i2 = 31.84°, 923 = 44.61°, ^13 ^ 7.82°. 
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